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Abstract: We consider the propagation of ultra high energy cosmic rays through 
{Nj ! turbulent magnetic fields and study the transition between the regimes of single 

and multiple images of point-like sources. The transition occurs at energies around 
E c ~ Z 41 EeV(B rms /5 /iG) (L/2 kpc) 3 ^ 2 ^/50 pc/L c , where L is the distance tra- 



versed by the CR's with electric charge Ze in the turbulent magnetic field of root 



mean square strength B rms and coherence length L c . We find that above 2E C only 
sources located in a fraction of a few % of the sky can reach large amplifications 
of its principal image or start developing multiple images. New images appear in 
pairs with huge magnifications, and they remain amplified over a significant range of 



energies. At decreasing energies the fraction of the sky in which sources can develop 
multiple images increases, reaching about 50% for E > E c /2. The magnification 
peaks become however increasingly narrower and for E < E c /3 their integrated ef- 
fect becomes less noticeable. If a uniform magnetic field component is also present it 
would further narrow down the peaks, shrinking the energy range in which they can 
be relevant. Below E ~ E c /10 some kind of scintillation regime is reached, where 
many demagnified images of a source are present but with overall total magnification 
of order unity. We also search for lensing signatures in the AGASA data studying 
two-dimensional correlations in angle and energy and find some interesting hints. 
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1. Introduction 

One of the fundamental open problems in physics is to understand the nature and 
origin of the ultra- high energy cosmic rays (UHECRs). The increased aperture of 
the next generation of CR detectors (e.g. the Auger Observatory) will be crucial to 
attempt to locate possible astrophysical point-like sources from the observation of 
clustering or anisotropics in the arrival directions. If the UHECRs are electrically 
charged, as would be the case in the likely situation that they are protons or heavier 
nuclei, a very important fact which has to be taken into account is their deflection 
by the magnetic fields present along their path. These deflections are expected to be 
quite large below the ankle (i.e. for E < 5 x 10 18 eV), and hence no CR 'astronomy' is 
possible below those energies since the information about the original source direction 
cannot be recovered from the observed arrival directions. On the other hand, since 
typical galactic magnetic fields (B ~ few /iG) are unable to confine CRs within the 
Galaxy for energies above the ankle, this supports the belief that UHECRs have an 
extragalactic origin, in agreement with the absence of strong anisotropies toward the 
Galactic disk being observed at those energies. The decreasing magnetic deflections 
for increasing energies lend us hope that a correlation among the arrival directions of 
different CR events and between those and the possible source directions will clearly 
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show up with increased statistics at energies E ^> 10 19 eV, when the CR trajectories 
progressively approach the rectilinear propagation regime. 

The astrophysical magnetic fields [|1|, 0, |3| may be classified under two broad 
categories: the regular and the random fields. Large scale regular fields may be 
produced by adiabatic compression of preexisting cosmic seed fields, by their ampli- 
fication by a dynamo mechanism or may be related to the action of a galactic wind. 
They are known to exist in our galaxy (and other spirals as well), where the field 
lines essentially follow the spiral structure, with reversals in direction taking place 
between neighboring arms. The existence of coherent fields in the halo or even larger 
scales is debatable, as is the presence of regular fields on cluster scales. The random 
component is present on galactic scales, probably originating from supernova explo- 
sions which feed power on typical scales of ~ 100 pc, which is latter transferred to 
smaller scales due to the turbulence present in the high Reynold's number interstellar 
medium. This leads to a Kolmogorov spectrum of magnetic field fluctuations extend- 
ing probably down to very small scales (~ 10 8 m and with an rms amplitude 
exceeding the typical values of the uniform field by a factor 2 3. In cluster's cores 
fiG fields have been observed, and their coherence length is believed to be at the kpc 
scale. For the intergalactic medium usually turbulent random fields are also assumed 
to exist, but with much larger coherence lengths (~ 1 Mpc). The amplitude of these 
fields is bounded from the non-observation of Faraday rotation of distant polarized 
radio sources, but it can be sizeable (~ 10 -8 -r- 10~ 7 G) if confined in thick sheets 
e.g. around the supergalactic plane. 

The existence of these magnetic fields will certainly affect the observable prop- 
erties of UHECRs (see e.g. [[5], || 0, [], |TU|). As we have shown in previous works 
[ p] , [121 , besides inducing a change in the arrival directions, the deflections caused 
by the magnetic fields can also lead to strong lensing phenomena, including the for- 
mation of multiple images and energy dependent magnifications or demagnifications 
of the CR fluxes. In those works we dealt with the effects resulting from the regu- 
lar galactic magnetic field, which were shown to lead to important consequences for 
E/Z < 30 EeV (where Ze is the assumed CR charge and 1 EeV = 10 18 eV). It is the 
purpose of the present work to extend the analysis of the magnetic lensing effects to 
the case of random magnetic fields, and consider also the possible interplay between 
the random and uniform components. 

One may distinguish between four different energy regimes when CRs traverse a 
distance L through a random field with coherence length L c . If 5 oc 1/E is the typical 
deflection suffered by the CRs, in the high energy limit, corresponding to 5 <C L c /L, 
CRs propagate almost rectilinearly and those arriving to an observer from a faraway 
source have all seen essentially the same magnetic field. When 5 ~ L c /L, we enter 
the regime in which multiple images of the same source can appear and CRs from 
different images have felt uncorrelated values of the magnetic field. We find that 
in this regime large magnification effects may be observable. For 5 > fewL c /L, a 
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regime with large quantities of secondary images appears, having the property that 
magnification effects tend to be averaged out, and essentially only a blurred image 
with characteristic angular size given by S will be observable. The fourth regime is 
that of spatial diffusion of the CRs, which is established when the gyroradius of the 
CRs is smaller than L c . For instance, for the galactic magnetic fields this happens 
below the ankle (for E/Z < 0.3 EeV), but may happen at much larger energies for 
UHECRs propagating across intergalactic magnetic fields. The energy dependence 
of the associated diffusion coefficients in this regime can lead to important changes 
in the slope of the CR spectrum, as is known to be the case for galactic CRs below 
the ankle. 

We will focus our work here in the systematic study of the first three regimes 
just described. Let us mention that focusing effects produced by random fields were 
previously noticed in the numerical study of ref. fL0"| , but were not analyzed in detail, 
while the first and third regimes were considered in refs. |7] in relation to the study 
of the time delays associated to bursting sources. 



2. The turbulent field 

We model the turbulent magnetic field as a Gaussian random field with zero mean 
and root mean square value B rms . This can be described by a superposition of Fourier 
modes as 

5*(x) = / ^^(k)e^ x+ ^) , (2.1) 

where the phases </>i(k) are random. If the turbulence is isotropic and homogeneous, 
the random Fourier modes satisfy the relation 

(^(k)^(k')) = ^^(27r) 6 5(k + k') , (2.2) 

where the projection tensor Pij = 5ij — kikj/k 2 guarantees that the field is solenoidal 
(V • B = 0) |Tj|. We will consider the general case of a power spectrum 

-L \*-'minl -'-'max) 

for 2ir/L max < k < 2tt / 'L min , and zero otherwise. This power spectrum is already 
normalized such that (|B(x)| 2 ) = B 2 . ms . The particularly interesting case of a Kol- 
mogorov spectrum (for which the energy density dE/dk oc k~ 5 ^ 3 ) corresponds to a 
spectral index n = 5/3. 

The correlation length of this field, L c , can be defined through 

/oo 
dL(B(0) • B(x(L))> = L c B 2 rms , (2.4) 
-oo 
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where the point x(L) is displaced with respect to the origin by a distance L along a 
fixed direction. The integral in the lhs of Eq. (|2.4j) can be computed using Eqs. (|2.1| ) 
and fl2.2|), and leads to 

T t B2{k) = LcBL ° (2 - 5) 

which can be used to express L c in terms of L m i n and L max as 

j 1 j- Tl 1 1 {J-i m i n I L m ax) (O R\ 

C ~ 7y max 1 77 77 \n-l l^""/ 
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J max i 



J max j 



In the case of either a very sharp (n ^> 1) or of a very narrow-band [Lmin ~ -^r 
spectrum, one gets L c ~ L max /2, while for a broad-band (L max 3> L m „) Kolmogorov 
(n = 5/3) spectrum, one gets instead L c ~ L max /5. 

The deflection in the velocity of a particle of charge Ze travelling a distance L 
through a magnetic field B, in the limit of small deflections, can be approximated 
by 

Ze 



S E 



L 

ds x B(V 



(2.7) 



For particles moving in a turbulent field the mean deflection vanishes. The relevant 
quantity is the root mean square value 



5 2 rms = (6 2 ) = Q-) 2 I'ds J Q L ds>(B ± (s)-B ± (s')} } (2.8) 

where B^(s) stands for the component of B in the direction orthogonal to the tra- 
jectory. This can be computed with the help of Eqs. ( p.l|) -( p73[) . In the limit of 
L ^> L max , corresponding to the distance travelled much larger than the maximum 
turbulence scale, the result can be written as 



y/2 E rmsV c E/Z 5 fiGy 2 kpcV 50 pc ' V ; 



where we have used Eq. ( [2.5|) to write it in terms of the coherence length. We have 
chosen to express its numerical value in terms of parameters relevant for the galactic 
magnetic field, which has a random component with rms strength of a few /iG and 
a maximum scale of turbulence of order 100 pc. 

Notice that the position at which the charged particle arrives forms an angle 
with respect to its initial velocity given by 



Ze 

which has a root mean square value 



[ L ds [ S ds'B ± (s') 
Jo Jo 



(2.10) 



meaning that after traversing a distance L in a turbulent field, the dispersion in the 
directions of the velocity vector is larger than that in the position vector as seen from 
the departure point by a factor y^- 
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3. The high energy regime 



As a beam of charged particles propagates through a magnetic field, its flux is fo- 
cused or defocused due to the differential deflections of neighboring paths, with the 
effect being larger at smaller energies. Eventually, the focusing effects may become 
sufficiently strong to produce multiple images of single CR sources. Heuristically, 
multiple images of a single source require CRs to reach from the source to their des- 
tination point traveling through uncorrelated values of the intervening magnetic field 
M. Consider a very distant source, so that its flux can be approximated as a beam 
of parallel CRs when they enter the region permeated by the magnetic field. The 
typical condition for the formation of multiple images is that two initially parallel 
paths separated by a distance of order L c may reach the same point after travers- 
ing a distance L in the magnetic field, i.e. that t]{L) ~ L C /2L, or equivalently 
SL/L C ~ y/3/2. Multiple images of a single CR source are thus a likely possibilty at 
energies E ~ -^E C) l with E c defined through 

$rms = ~f^~j~ ■ (3-1) 
Hz L 

Its numerical value reads 

Notice that the heuristic argument above applies if the magnetic field has one 
dominant length scale only. If the turbulent field has a broad-band spectrum, this 
heuristic argument should be applied to the effect of each individual wavelength bin. 
Small wavelengths can actually lead to multiple image formation at higher energies 
than long ones if the amplitudes of the Fourier modes of the magnetic field do not 
decrease too fast at small scales, as will be discussed further below. 

A detailed study of the focusing effects of the magnetic field can be performed 
following the trajectories of neighboring paths in a beam. The amplification of the 
flux of a CR source, given by the relative change in the area of the cross section 
of an initially parallel beam, can be written borrowing the formalism familiar from 
gravitational lensing [|14j] as 

A = — — — . (3.3) 
|(1 - k) 2 -Y\ 

Here the convergence k describes the isotropic (de)focusing effect while 7 is the shear 
that distorts the shape of the cross section of the beam. It is useful to decompose 

x If the source rather than being practically at infinity were instead at a distance L within the 
spatial extension of the magnetic field, the typical focusing condition is that at the midpoint the 
separation between alternative paths be of order L Cl i.e. r)(L/2) ~ L c /L, which corresponds to 
E ~ E c j\f%. This is 2\/2 times smaller than in the case of a very distant source. 
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the amplification along a particular direction of observation as A = A1A2, with 
A\ = 1/(1 — k — 7) and A2 = 1/(1 — k + 7) measuring the relative stretching of 
the beam along the so-called shear principal axes. This means that a beam with an 
initially circular cross section with diameter L develops an elliptical cross section, 
with one of the principal axes having its length changed by ALi = Lq(—h — 7) and 
the other one by AL 2 = L (— K + 7). As the deflections due to the random field have 
no preferred direction 2 , this means that (ALi) = (AL 2 ) = and (AL 2 ) = (AL 2 ). 
Since ^(7) ~ — [ALi + (— )AL 2 ]/2L , it is then easy to see that (k) = (7) = and 
that (k 2 ) = (7 2 ). 

In the limit of small deflections the isotropic focusing of an initially parallel beam 
of CRs is given by 

1 Ze r L 

k = -— (L-s)VxB.ds. (3.4) 

2 E Jo 

Using eqs. fl2.1|)-(|2~3T) we can calculate (k 2 ) = K 2 ms , with the result 

_/-_L^R L3/2 JrmsL ,E C 
Krms — S jj, D rms ^jj- ~ S> — E ' 

where we defined the numerical coefficient ( such that 

2 - n A (n-l) 2 l-(WW n nfi , 
) n{2-n){l-{L mm /L max y-^ ' 1 ' } 

In the case of either a very sharp or of a very narrow-band spectrum, one gets 
C = tt/v12 ~ 0.9. For a broad-band Kolmogorov spectrum, one gets instead 
( ps 0.8 [Lmaxl L-min) 1 ^ ■ Notice that, contrary to the case of deflections, small 
wavelengths have a significant effect upon focusing of charged particles for spectral 
indices n < 2. Indeed, if n < 2 the energy at which n rms becomes of order unity, 
and thus strong lensing effects become very likely, increases by a factor of order 
(Lmax I Lmin) 1 ~ n ^ 2 compared to the case of a very sharp spectrum. 

In the high energy regime, i.e. for E > ( E c , k and 7 are small and they can 
be approximated by uncorrelated random Gaussian variables with zero mean (see 




Eq. (|3.4|) and take into account the fact that B is Gaussian) and dispersion K rms . 
The probability distribution of the amplification of the beam in different observing 
directions, assuming that the observer is located at the center of a sphere of radius 
L filled with a random field of constant strength B rms is then given by 

d 2 p 



d/td7 



G K (0, K rms )Gy(0, K rms ), (3-7) 



2 Contrary to the gravitational lensing case, for magnetic lensing the convergence can also be 
negative, and the shear can have any sign since we are not making any distinction between the two 
principal axes. 
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with G x (x, a) = exp[(x — x) 2 /2o 2 }/y2iia. 

More interesting is to study the probability p s that the flux from a distant source 
be seen by the observer with a given magnification. As long as we can neglect the 
presence of multiple images, which is a reasonable approximation for E > 2E C , this 
one is related to the probability in Eq. ( |3.7| ) by a factor 1/A, since the magnification 
is just the ratio between the observed solid angle to that actually subtended by the 
source (because the surface brightness is conserved). Hence, one has that in the high 
energy domain 

d 2 p s d 2 p 



:i--) 2 -7 2 i 



dud'y dud'j 1 v ' ^ 

One can then compute the probability that a source be magnified above a certain 
threshold A as 



P S (A > Ac) ^ J d 7 J d K 



d 2 p s 
dudj 



This gives that e.g. for E 
by a factor larger than A 
of five. 



(3.9) 

(LMl . 

3E C around 6% of the source directions are magnified 
2, while only 1% are magnified by more than a factor 



0.25 



I 0.15 




0.05 



Figure 1: Fraction of the source directions for which multiple images would appear vs. 
E/E c . 



The flux from a CR source diverges after traversing a distance L in the magnetic 
field for those directions such that (1 — k — 7)(1 — k + 7) = 0. In that case, the 
source is located on top of a caustic of the magnetic field configuration relative to 
the specific point of observation. Multiple images of a given source are visible at 
energies below that at which a caustic crosses the source location. The additional 
images exist in pairs of opposite parity, and for the inverted image one has that 
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(1 — /«) 2 — 7 2 is negative. We can then compute the probability that a given source 
has multiple images as 

oo 1+7 ^ 2 

Pmuit = P a (\l ~k\< | 7 |) = 2|d 7 J dK^p (3.10) 

1-7 

It has to be noticed that this expression is reasonably accurate as long as the fraction 
of sources with more than three images is negligible (which corresponds to approxi- 
mately Pmuit < 25%). The result is plotted in Figure [I], and we see that for E ~ 2E C 
around 3% of the sources will have multiple images and this fraction rises to ~ 20% 
at E ~ 1.25E C . 




E/E c 



Figure 2: Mean number of images of a source vs. E/E c . 

For decreasing energies, the number of images associated to a source increases. 
As will be shown later on in Section ^ (see Eq. ( |6.3|) ), this number can be computed 
as the mean of A~ l over all the observing directions 9. Taking the mean over many 
independent directions is actually equivalent to take the mean over realizations of 
the random field, i.e. over different values of k and 7. We then obtain 

{n) = h I d2M_1 = I dKd7 o^ l(1 _ K)2 ~ l2y (3 - n) 

This expression is plotted in Figure as a function of E/E c . We see that (n) 
takes values 1.06, 1.8 and 4.8 for E/E c = 2, 1 and 1/2, showing that multiple imaging 
is of paramount importance for E < 2 E c . 

4. Numerical results 

We have implemented a numerical code to propagate charged particles within a real- 
ization of a random magnetic field, obtained as the superposition of iV independent 
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waves, with random directions of the vector k and direction of B(k) randomly chosen 
in the plane orthogonal to k, so that the condition V • B = is automatically fulfilled. 
The modulus of each mode |B(k)| is drawn from a Gaussian distribution with zero 
mean and dispersion given by Eq. ( |2.3|) with n = 5/3 (Kolmogorov spectrum). The 
code allows to calculate not just deflections of trajectories but also flux amplifica- 
tion along them, by following two additional nearby particles, using the formalism 
developed in ref. In practice, CR trajectories that reach the Earth are found 

backtracking trajectories of particles with opposite charge from the Earth out to a 
sphere of radius L. 



The code results agree with the analytic expressions in Eqs. (|2.9| ) and (|3.5| ) for 
the rms deflections and amplifications in the small deflection limit for a wide range 
of parameters. 

A useful way to visualise the lensing properties of a magnetic field is to plot, 
for a regular grid of CR arrival directions on Earth, the incoming directions they 
had as they entered the region permeated by the field. These 'sky sheets' [^TJ are 
quite smooth, regular and single-valued at energies E ^> £E C , meaning that only 
one image of each source is visible at high energies. Its flux will be demagnified or 
magnified in proportion to the amount of stretching or compression of the sheet in 
the position that corresponds to the direction with which CRs enter the magnetic 
field. For decreasing energies the sheet becomes eventually significantly folded due 
to increased relative deflections. CRs that enter the magnetic field from a direction 
where the sheet is folded are seen on Earth as coming from all those directions in the 
grid that overlap in that point of the picture. The additional images of the source 
are visible at energies below the one at which the fold (which moves as as function 
of energy) crosses the source position. 

We display these effects in Figure ||, for three representative values of E/Z. The 
figure corresponds to the case of a field strength B rms = 5 /iG, a sharp spectrum 
(L max = L min ) with turbulence scale L max = 100 pc (L c = 50 pc), and the distance 
traveled by the CRs within the field is L = 2 kpc. These values are representative of 
the random component of the galactic magnetic field, and we shall use them along 
the paper to illustrate numerical results. The effects are rather universal for different 
sets of parameters in terms of E/E c , but the angular scales involved also depend on 
L/L c . The simulations were performed with N = 300 independent waves to generate 
the random field, and the step with which the trajectories were numerically followed 
was 1/5 of L c . 

It is clear from the left panel in Figure § that folds already cover a small but non- 
negligible fraction of the sky, of order a few %, at energies around 2E C . This fraction 
grows considerably at energies around E c (center), in agreement with the values 
obtained for P mu it in the previous Section. At energies around E c /2 (right panel) 
folds on top of folds have already developed, which imply sources with more than 
a pair of additional images. The folds cover about half the sky at energies around 
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E/Z = 82 EeV = 2 EJZ E/Z = 41 EeV = EjZ E/Z = 20.5 EeV = EJ2Z 




Figure 3: 'Sky sheets': projection of a regular grid of CR arrival directions on Earth 
onto the direction they had when they entered the region permeated by the magnetic field. 
Sources located in regions where the sheet is stretched are demagnified, and viceversa. 
Sources located in regions where the sheet is folded have multiple images. The parameters 
are B rms = 5 fj,G, L max = L min = 100 pc, L = 2 kpc (E c = 41 EeV). 

E c /2, and almost the full sky at energies between E c /3 and E c /4, meaning that 
whatever source one choses will have already developed multiple images at energies 
larger than E c /4. 

In Figure |] we illustrate the effects of small wavelengths in the turbulent random 
field through the analogous sky sheets for the case in which the rms strength of the 
field is the same as above (B rms = 5 /iG), the distance traversed is also L = 2 kpc, 
but the turbulence scales extend from L max = 100 pc down to L min = L max /10 with 
a Kolmogorov spectrum. In this case L c « 25 pc and E c « 58 EeV. We display the 
results for energies twice as large as in Figure ^ corresponding to an enhancement 
above E c by a factor {L max / L min ) 1 ^ . We used in this case N = 3000 independent 
waves to realize numerically a sufficiently random field. 

The comparison between Figures |3] and [| makes evident that the short wave- 
length modes of the turbulence lead to folds at larger energies than the long wave- 
length modes, but their angular scale is smaller (notice that the linear scales in 
Figure ^ are three times larger than in Figure |3|). However, as we shall discuss later, 
the integrated magnification effects upon a source depend not only on its proximity 
to a fold but also on how fast the fold moves across the source location, which is 
basically determined by the long wavelength modes, and this may average out the 
effects associated to the short wavelengths. 

The numerical code can also follow the change in apparent position and the flux 
amplification of a given source image as a function of energy, exploring at each energy 
step a neighborhood of arrival directions from the previous step, and backtracking 
it until it reaches the appropiate direction to the source with the required accuracy. 
The arrival directions of multiple images at some fixed energies found in the data 
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E/Z = 164 EeV = 2 (L max / L mi J' /e EjZ E/Z = 82 EeV = (L ma jL mi J' /e EjZ 



■5.8 I J .5.8 




Figure 4: Sky sheets with the same parameters for the magnetic field as in Figure |3| but 
with a Kolmoe orov spectrum for the turbulence, extending from Lrnax 

= 100 pc down to 

Lmin = L max /10. The energies are two times larger than in the left and center panels 
in Figure ||. Notice that the scales are also different, to better visualize the existence of 
structure on smaller angular scales. 



used to plot the sheets above are used as starting points to follow them in energy. 
In Figure |^ we display the amplification as a function of energy for two illustrative 
diverse behaviours. 

The solid line in Figure [5] is illustrative of the effects that sources located in only 
a small fraction of the sky, of order 2% or less, undergo due to magnetic lensing. 
This source is close to the location where a fold in the lens mapping develops at 
relatively high energies. Its principal image is highly magnified over a significantly 
broad energy range. Multiple images are visible at relatively high energies, with 
magnification peaks that have strong integrated effects. Indeeed, the magnification 
of the secondary images behaves around the energy Ef at which the fold crosses the 



source location as A = A E /yl — E/Ef |i~2"| , with a coefficient Ae that we found to 
be typically of order but larger than unity if Ef is above 2 E c . Since the magnification 
integrated over an energy bin of order 10% of the peak location is around 12 Ae (see 
next Section), we conclude that quite strong features are certainly imprinted upon 
the observed spectra of these sources, albeit with a small probability. The principal 
and the first pair of secondary images are demagnified at decreasing energies, but 
many additional images (not depicted here) will eventually become visible, although 
with much narrower peaks. 

The dashed lines in Figure |5] are illustrative of a more likely situation. The first 
pair of secondary images of this source develops at an energy around E c . At relatively 
high energies the source location is away from a fold, in a stretched region of the 
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40 60 80 100 120 140 160 180 200 

E/Z (EeV) 

Figure 5: Amplification vs. energy for the principal and first pair of multiple images of 
a source location close to where a fold develops at energy above 2 E c (solid lines) and for 
a source location which is first crossed by a fold at energy around E c (dashed lines). The 
parameters are the same as in Figure || (E c = 41 EeV). 

sky sheet, and thus its principal image is demagnified. The magnification peaks of 
the first pair of secondary images are relatively narrow, but still lead to significant 
integrated effects, since typically Ae is around 0.5 if Ef is around E c . We find that 
typical values of Ae remain above 0.1 (below which no significant effects would be 
observable unless the experimental energy resolution were better than 10%) down to 
energies Ef around E c /3. 

At energies below E c /3 it becomes very likely that sources develop a high num- 
ber ((n) > 10) of secondary images. It is thus unpractical to trace each individual 
image apparent location and amplification as a function of energy. A ray shooting 
technique becomes instead more appropriate. In such simulations, a large number of 
antiparticles are thrown isotropically from Earth, and those that after traversing a 
distance L point to a direction closer than a required accuracy from that to a fixed 
'source' are recorded. The ratio of this number to the one that would have been 
obtained in the absence of magnetic deflections is just the corresponding magnifica- 
tion of the total source flux, summed over all images. For a fixed number of rays 
shooted, a smaller source size degrades the precision with which demagnifications are 
recorded. A large source size instead prevents recovering large magnification peaks, 
due to averaging effects. 

In Figure || we plot the result of evaluating through such ray shooting technique 
the amplification of 1080 consecutive source locations of angular diameter 1/3 of 
a degree each. The parameters are the same as in previous Figures. We plot the 
magnification of each image sorted from higher to lower values, with the horizon- 
tal axis rescaled to 1. Thus, the value of the abscisa for a given amplification A 
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E/Z=250 EeV 
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0.2 0.4 0.6 0.8 1 

f 

Figure 6: Fraction of sources (in the horizontal axis) with amplification higher than a 
given threshold (in the vertical axis), obtained from the ray-shooting simulation described 
in the text, for representative values of E/Z. The fit to Eq. ( |3.9[ ) is also shown for the 
highest energy. The magnetic field parameters are as in Figure |3| 

corresponds to the fraction of sources that have amplification larger than A . The 
results are displayed for three representative values of E/Z. At energies E ^> E c 
about half the sources are sligthly magnified while the other half are slightly de- 
magnified, in comparable proportions. The probabilities here are in agreement with 
Eq. (|3.9|) , which is also plotted for comparison. At decreasing energies more sources 
are demagnified rather than magnified, and there is a significant enhancement of the 
large magnification tail. The large magnification tail reaches its maximum strength 
around E ~ Z 50 EeV ~ E c . The fraction of sources demagnified also reaches its 
maximum, of order 2/3, at E ~ E c . At energies below E c the curve starts to level 
off again, as we enter some kind of scintillation regime where each source has a very 
large number of demagnified images with total amplification of order unity. Only 
the very few sources that have images with extremely narrow peaks precisely at the 
energy under consideration make a small contribution to the high magnification tail 
(but which is here somewhat suppressed by finite source size effects). 

In the next two Sections we show how these issues can be addressed analytically. 

5. The appearance of new images 

As we have seen from the previous simulations, new images of CR sources can appear 
when the mapping from the source's plane to the observer's plane becomes multiply 
valued or, in other words, when the observer's sky becomes folded when projected 
into the source's sky (Figures [3] and f|). This pictorial view of a folded sky is quite 
useful to study the general properties of the images in the regime of multiple imag- 
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ing. To see this let us consider the deflection of CR trajectories after traveling a 
distance L through a turbulent magnetic field and let us analyse the properties of 
the mapping between the incident directions (/?) and the observed ones (0). At large 
energies, trajectories are straight lines, and the mapping is the identity, but the typ- 
ical deflections increase with decreasing energies and make this mapping non-trivial. 
Since directions separated by more than Lj L probe uncorrelated values of the mag- 
netic field, this means that they suffer uncorrelated deflections. These uncorrelated 
deflections have a random distribution of directions and a typical amplitude given 
by S rms . The two-dimensional deflection process can be thought as the superposition 
of one-dimensional deflections in two orthogonal directions, and we show in Figure |7| 
a picture of a plausible mapping for one of the directions. 




L c /L 

Figure 7: Pictorial view of the mapping between the arrival direction of cosmic rays at 
the observer (6) and the incident direction from the source (/?). 

The general properties of this mapping can be simply understood considering a 
network of directions separated by L c / L and assuming that these points are deflected 
with a certain (energy dependent) amplitude in either direction. Two neighboring 
points in this network can then be deflected towards the same or in opposite direc- 
tions. If they are deflected towards the same direction, essentially only the position 
of the image changes. On the other hand, if two neighboring directions are deflected 
toward each other (so that A/3 < AO), a caustic will eventually form in between 
them (when d/3/dO = 0) as the energy decreases, leading to large amplifications of 
the flux and to the formation of multiple images. Since in this picture neighboring 
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directions have a probability 1/4 of approaching each other, the typical separation 
between the directions where caustics will form is Af3 = 4L C /L, in agreement with 
the results of the detailed simulations. When neighboring points are deflected in 
opposite diverging directions (so that A/5 > A9), a low magnification region results. 
This low magnification regions will appear always in between the caustics. 

Around a direction where a caustic forms, which is the most interesting case due 
to the possibility of having large magnification effects, the mapping in the direction 
orthogonal to the fold can be approximately described by a cubic polynomial. Defin- 
ing y = {j3 — P )L/L C and x = (9 — 9 )L/L C (where /3 is the direction where the 
caustic forms and 9q its image, i.e. the corresponding critical line in the observer's 
plane), the mapping can then be written as 3 

y = ax + bx 3 . (5.1) 




Figure 8: Mapping near a caustic for three different values of E*/E = 0.5 (dashed), 1 
(solid) and 1.5 (dotted). 



The coefficients a and b are functions of the deflections 6(E). For E oo, we 
have a — > 1 and 6^0, and it is natural then to take that a — 1 oc 5 oc 5 oc E~ x . 
The caustic forms at x = when dy/dx = 0. This means that a can be written 

3 For simplicity we omit a quadratic term, what corresponds to the assumption that the positive 
and negative deflections generating the folds have similar magnitudes. The generalisation to the 
asymmetric case can easily be done by assuming that the symmetry point of the fold in the source 
plane, /?o, moves with energy. The leading correction would then be obtained replacing [3q — > 
(3o + D(E — E*), with D a constant. The main effect of this change is to modify the energy gap 
between the first bump in the spectrum and the appearance of the first fold peaks, as well as the 
energy width of these last (see below). 
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as a = 1 — E*/E, where E* denotes the energy at which the caustic forms, and 
6 = CE*/E with C being a constant. The value of C can be inferred from the 
fact that the slope dy/dx has to be unity at some point near x ~ ±1/2. Hence, 
we should have that typically C ~ 4/3. The mapping y(x) is depicted in Figure |] 
for different values of the energy, showing the formation of caustics as the energy 
decreases below E*. Notice that the energy E* will be different for different source 
locations, since it is determined by the energy of formation of the fold nearest to 
the source chosen. An average value of E* can be estimated from the fact that 
y(x = 1/2) = 1/2— E#/3E ~ 1/2— 5L C /L, and as the rms value of the one dimensional 
deflection 5 is just 5 rTns /\/2 this leads to (E*) ~ (3/y/2)E c . 

For E < E*, a couple of folds appear that move appart as the energy decreases. 
Their location in the source plane, ±y/, can easily be found from the condition 
dy/dx\ Vf = 0, leading to 

IE* ( E \ 3/2 , . 

^M'-k) ' (5 ' 2) 

Particles from an incident direction y s will be observed from the direction(s) x 
that solve the cubic equation (|5.1|). The solution is unique (just one image of the 
source) for E > Ef, where Ef (which is smaller than E*) is the energy for which the 
fold location coincides with the source direction. This means that Ef is just obtained 
from the relation 



Vs 3 E f V E. 
The location of the image in this case is given by 

'• :; i -q+ Jo) 1 + (~q 



3/2 



1/3 



(5.3) 



(5.4) 



p 3 + q 2 . 



where we introduced q = —y/(2b), p = a /(3b) and Q 

The amplification of an image (in the direction orthogonal to the caustic) is given 

by 

d9_ _ fdy 
d/3 \dx / 

Using Eq. (|5.4|) we then get 
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(a + 36a; 
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(5.5) 
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(5.6) 
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Notice that in the demagnification regions, the deflections will just have the opposite 
sign than what was assumed before, and hence a similar expression will hold for the 
amplification but changing E — > — E in the expressions for a and b. 

When the source is located within the folded region, i.e. for \y s \ < yf, there are 
three different images, whose positions can be written as 

a 



x 



(k) 
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with k = 0,1,2. Here k = corresponds to the principal image (the one already 
present at high energies) , while k = 1 , 2 are the pair of images created when the fold 
crosses the source location (i.e. for E < Ef). The amplification of these images is 
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3y s E/E* 
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In general the total amplification will be the product of the amplification in the 
direction transverse to the fold, A±, times the amplification in the direction along 
the fold, An (that can be taken as constant in first approximation). Then, one has 



that A to t = A\\ J2k l-^j^l- We notice that a negative value of A±' would just mean 
that the corresponding image is inverted, i.e. that it has negative parity. 



(k) 




Figure 9: Amplification of a point source for different values of y s 
1/80 (long dashes), 1/40 (solid) and 1/4 (dotted). 



1/160 (short dashes), 



In Figure |9] we show the transverse magnification of the images of a source for 
different values of y s (i.e. for different values of the distance from the source to the 
symmetry point of the fold in units of L c /L). When plotted in terms of E/E*, these 
curves only depend on the value of y s , with smaller y s leading to more pronounced 
magnification effects. Notice the similarity of Figure ^| with the energy dependence 
of the amplification in the numerical examples in Figure [|. 

Let us now first concentrate in the magnification of the principal image of the 
source. We find that for small values of y s (y s < 0.1), sizeable peaks are observed. 
The theoretical amplification obtained from Eq. (|5.6|) provides an excellent fit to the 
amplification of these peaks obtained in numerical simulations, as is exemplified in 
Figure [TO. Notice that the value of A max determines y s , while the location of the 



peak fixes E*. 

The relation between the peak magnification achieved, A max , and the source 
location y s is depicted in Figure [IT] (solid line). This curve is accurately fitted, for 
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Figure 10: Amplification of the principal image of a point source computed numerically 
(+) and fit using the theoretical model. 

y s < 0.1, by the expression 



as is also illustrated in that Figure (dashed line). 

The peaks in the magnification of the principal image for y s < 0.1 are quite 
wide, having typical widths of 10 -j- 30% of the peak energy. This is very interesting 
because it means that these peaks have a width comparable to the energy resolution 
of UHECR experiments, and should then be in principle resolvable. 

We can also estimate the probability that the peak magnification of the principal 
image be larger than a given value A$. This equals the probability that the source 
be at a distance to the caustic formation line smaller than y Q L c /L, with y related to 
A through Eq. (|5.9|) . Assuming that a rectangular network of caustics is formed in 
the sky separated among them by 4L C /L (i.e. by Ay ~ 4), this probability will just 
be equal to the fraction of the sky covered by a network of strips of width 2y L c /L, 
and separated among them by 4L C /L, which is just equal to yo, i.e. 



with the last equality being valid for Aq > 2. This expression was checked numeri- 
cally by following the principal images of a thousand randomly located sources and 
computing their amplification, and it indeed agrees very well with the numerical 
results obtained. These results imply that there is for instance a one percent proba- 
bility for the principal image to be magnified by more than a factor eight, while the 
probability for it to be magnified by more than a factor of three is ~ 5%. 

Turning now to discuss the magnification of the new pair of images appearing 
below Ef, we see from Figure |^ that they have associated divergent peaks, which for 
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Figure 11: Maximum amplification of the principal image of a point source as a function 
of y s (solid line). The dashed line corresponds to the fit described in the text. 



increasing values of y s appear at smaller values of Ef/E* and become increasingly 
narrower. We can obtain an analytic approximate expression for the amplification 
of these peaks making a Taylor expansion of Eq. (|5.2j) around the fold location 
Xf — ±J—a/3b. This gives (for the fold at positive values of y) 



V-Vf+ oVf( x ~ x f) 2 + ■ 



with 



y.f 



d 2 y 



dx 2 



E* E 
~E\ ~ EZ' 



(5.11) 



(5.12) 



Deriving the expression Eq. ([5.11 ) we obtain the amplification of the flux of CRs of 
energy E coming from a source in the direction y s , which is 



A - = {fx) ^±(2«J(y.-y/) 

V / u.. 



-1/2 



(5.13) 



Furthermore, since the fold moves with energy according to Eq. ( |5.2[ ) while the source 
is at a fixed position y s , we can write yf ~ y s + (dyf /dE) Ef (E — Ef), with 
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Combining these expressions, one gets for energies close to Ef that 4 

A ■ Ae 

A i ~ : 



l-E/E/ 



(5.14) 



(5.15) 



4 Notice that in this approximation the magnification of the two images is similar and their 



difference only shows up if we include additional terms in the Taylor expansion in Eq. (5.11) 
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showing the characteristic 1/yl — E/Ej divergence and with the coefficient Ae being 



A, 



2y f Ef dE 



-1/2 




- 1 1 + 2 



-1/2 



(5.16) 



It is important to notice that the coefficient Ae depends on how fast the magnification 
changes with the observation direction (y"f~ l ) and also on how fast the fold position 
moves with energy (given by dyf/dE). 

To find out what are the possible observable effects of these magnification peaks, 
it is useful to introduce the amplification integrated in an energy bin around Ef (e.g. 
between 0.9Ef and Ef), since this would be indicative of the potentially observable 
signals in a realistic experiment once its finite energy resolution is taken into account. 
This has the further advantage that the integrated magnification becomes finite. 
Hence, we define 



Ai n t = 
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[ Ef dE J2 12-6^1 
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Figure 12: Integrated amplification Ai n t in an energy bin of 10 % of Ef as a function of 
the source position y s (solid line). The dashed lines correspond to the fits described in the 
text. 



In Figure [L2] we plot the integrated amplification Ai nt as a function of y s . The 
decrease in A int for increasing values of y s is mainly due to the fact that the peaks 
become increasingly narrower for smaller Ef/E*. The long dashed curve corresponds 
to A int = 3.3y~ 38 , and provides a reasonable fit for y s < 0.2, while the short dashed 
curve corresponds to A int = 1.9yJ 07 , fitting the results for 0.2 < y s < 3. From 
this we can obtain the probability that the integrated amplification Ai nt exceeds a 
given value A , which equals the probability that y s be smaller than the value yo 
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obtained from the condition that Ai nt = Aq. Using the previous fits we then get (for 
AE = O.lEf) 

P(A mt > Ao) = P(y s < y ) = yo * | 2 .5Ao 143 A < 6 (5 ' 18) 

This results imply that there is for instance a probability of ~ 50% that the secondary 
images lead to an integrated amplification A int larger than three, while A int > 7 with 
~ 10% probability. Hence, the effect of these peaks can be in principle quite relevant. 

It is important at this point to consider what would be the implications of having 
in addition to the random magnetic fields also a regular one, coherent over a scale 
L reg 3> L c . If both components have comparable strengths, as happens in the Galaxy, 
magnification effects associated to the random one will manifest at higher energies 
than those of the regular one. This is because the lensing depends essentially on 
field gradients, which are enhanced when the field variations occur on smaller scales. 
However, the possible lensing signatures produced by the random field will be affected 
by the presence of the regular one. The most important effect will be related to the 
way in which the regular field changes the "motion" of the folds with decreasing 
energies, i.e. the factor dyf/dE in Eq. ( [5.16 ). The motion of the fold is essentially 



due to the change in the typical deflections with energy, and hence for purely random 

fields it is given by \dyf/dE\ ~ |d<5 / ~ 5/E. On the other hand, in the presence 

of a regular field, the typical deflections will be given by 

B± L reg 10 EeV 
Ag~0.1^ - reg — — , 5.19 
/iGkpc E/Z v ' 

with B± being the typical strength of the regular field in the direction orthogonal 

to the CR trajectory. From Eq. (|5.16|) we then see that the magnification peaks 

associated to new image pairs should become narrower, with their integrated effect 

being suppressed by the factor 




\/^( a65 if^^l • (5 - 20) 

where L > L reg is the distance traversed across the random field. For typical galactic 
parameters (i.e. B± ~ B rms /2, L reg ~ 3 kpc and L c ~ 100 pc), this gives 

J— -0.28 — . 5.21 

V A8 \3 kpc/ v ; 

Let us finally mention that something similar happens when we look at the peaks 
associated to the short wavelength modes of the turbulence. These modes can in 
principle produce peaks on very small angular scales (and even at somewhat larger 
energies that the long wavelength modes), but the motion of the small scale fold 
will be determined by the long wavelength modes (or the regular field if that one 
is present). This will make the small scale peaks extremelly narrow and hence less 
noticeable. 
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6. The scintillation regime 



As we have seen, the appearance of folds in the mapping from the image plane to 
the source plane is associated to the formation of pairs of new images. Since for 
decreasing energies the folds cover an increasing fraction of the source sky and also 
new folds are continuously generated, this means that the average number of images 
of a source increases steadily for diminishing energies. The transition between the 
regimes where only one or a few images exist and that in which many images are 
present may be modelled with the following simplified picture, which captures the 
main processes involved in the multiple image formation. Let us assume that first 
at a single energy E% a two dimensional rectangular grid of caustics forms, with 
separation among them A/3 ~ AL C /L. The energy E\ would just be the mean energy 
of appearance of the folds, i.e. (E+). For E < Ei, the two folds in each caustic 
become separated among them by a typical angular distance 2yfL c /L, with yf given 
by Eq. fl5.2| ) with E* — > E\. Since there are two additional images (a total of three) 
in the regions covered by a fold, while at the fold intersections there are 3x3 = 9 
images, it is then easy to show that the average number of images for randomly 
located sources is 



When the folds become sufficiently wide so that in each face of a fold there are 
directions probing uncorrelated magnetic field values, a second generation of folds 
can then be generated. This can be modelled by means of a new network of caustics 
appearing at an energy E 2 (corresponding typically to the energy at which yf ~ 1). 



This process should then repeat itself again at lower energies, leading to an expo- 
nential growth in the mean number of images of the CR sources. As an example, 
Figure [13] shows the numerical results for (n) and the fits obtained with the previous 
expression corresponding to one (long dashes) or two (short dashes) generations of 
folds appearing at energies E\ = 85 EeV (approximately 2 E c ) and E 2 = 20 EeV, 
with the results being indeed quite satisfactory. 

The numerical result for (n) in Figure |13| was obtained using the property that 
the amplification satisfies A(6i,6 2 ) = J~ l , where J is the Jacobian of the mapping 
between the observer's (#j) and source's (/%) coordinates. Hence, one has that 



i.e. that the average of the inverse magnification in the observer's plane is just the 
average of the image number in the source plane, and the first can be obtained 




(6.1) 
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Figure 13: Average number of images of a source vs. E, and theoretical fits (described 
in the text). In this example E c = 41 EeV. 



computing the magnifications for a dense grid of directions isotropically distributed 
around the observer. The theoretical expression obtained in Eq. (|3.11|) reproduces 
these results to better than 20% for E > E c /5, below which the appearance of 
multiple folds makes the assumption of Gaussian distributions for k and 7 certainly 
no longer valid. 

We can also use this simplified picture of a first network of caustics forming at an 
energy E\ to estimate the fraction of the sky in which sources have multiple images, 
as a function of energy. With arguments similar to those that lead to Eq. ( p.lp , we 
estimate f s ~ yj — y|/4. Using E\ = 85 EeV we get that more than 20% of the sky 
has multiple images around E = E c , and practically all the sky has multiple images 
already at energies between E c /3 and E c /4, as already found before by different 
means. 

Figure |l~3| shows that already for E = E c /3 there is an average number of images 
(n) ~ 8, and that for E = E c /5 this number has increased to (n) ~ 30. This large 
number of images is due to the continuous creation of image pairs, which appear 
largely magnified but the width of the peaks become increasingly narrow as the 
energy diminishes, and hence their integrated effect is reduced. As a result one 
reaches a regime having many demagnified (i.e. with Ai nt < 1) images of every 
source. These images will be spread over a typical angular scale which is 

not necessarily large if L c -C L. Furthermore, the total magnification of the many 
images (averaged in energy bins) will become of order unity in this regime. This can 
be understood from the property that (we are considering here an observer at the 
center of a spherical region filled with random fields of constant B rms ) 
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If one is in the regime in which all source directions have associated a large number of 
images, one may consider that all points in the source plane are essentially equivalent, 
and hence this leads to (Y,images A4) — 1- 

The properties of this 'scintillating' regime, which in some respects is reminiscent 
of the twinkling of the stars produced by the atmospheric turbulence, are illustrated 
in Figures [14] and [3j|, where we present the results of a (cosmic)ray-shooting. In 
this simulation a large number of anti-particles were thrown isotropically from the 
'detector', and those that after traversing a distance L point to a direction closer 
than 1/3 of a degree from that to a fixed 'source' were recorded. The ratio of 
this number to the one that would have been obtained in the absence of magnetic 
deflections is just the corresponding magnification. This was repeated for different 
energies and the results are plotted in Figure [TJ]. Superimposed in the same figure 
are the magnifications of the principal image of the source (the one visible at the 
highest energies) and of the first few pairs of secondary images. These were evaluated 
numerically with the method described in Section f|, which consists in tracking the 
trajectories of three nearby particles for each image. 




E/Z (EeV) 



Figure 14: Magnification vs. E/Z of an extended source from a ray-shooting simulation 
(solid line). Also shown are the magnifications of individual images (dashed lines) of a 
point source in the same location. The parameters of the simulation are the same as in 
Figs. ^, H and o, corresponding to E c = 41 EeV. 



Notice that in the ray-shooting simulation the divergences in the peaks are 
smoothed by the finite size of the source. A similar smoothing would also result 
for a point-like source due to the finite energy resolution of realistic detectors. The 
peaks associated to the first few pairs of images are clearly noticeable. As the en- 
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ergy decreases, the peaks become increasingly suppressed in width, and there is a 
progressive transition to the regime with (X^i) — 1- 

Notice that in this simulation E c ~ 41 EeV, which is about the energy at which 
the first peak is located. The source in this example is thus a rather generic one. 
Sources located in about 20% of the sky should display peaks similar or stronger 
than these. The impact of the peaks is generally enhanced by the fact that the 
fluxes are usually demagnified above the peak energy. Sources in the rest of the sky 
would display somewhat narrower peaks. On the other hand, a smaller fraction of 
source locations would have magnification peaks at higher energies associated to the 
principal image, instead of having it demagnified as in the example of Figure [14. 



In Figure [15] we show how the different images will look like at different energies. 
Each point in the top panel of Figure [15] represents 1/10 of the unlensed flux of 
the source, while in the bottom panels each dot represents 1/100 of it. Before the 
appearance of the first pair of images the principal one is only slightly displaced (it 
is assumed to be at the origin at very high energies) and demagnified (the source 
is not very close to the fold formation location). At lower energies, new magnified 
images appear and they are spread over the angular scale S rms oc 1/E. At the 
smallest energies considered, there are very many demagnified images but with total 
magnification close to unity. 

Notice that 5 rms becomes of order unity, and then multiple images of a source 
cover a significant fraction of the sky, at energies around and below E c /N c (see 
Eq. fl3-ip ), where N c = L/L c is the number of incoherent domains of the magnetic 
field traversed in a straigth line (in our example E c = 41 EeV and N c = 40). Spa- 
tial diffusion sets in at somewhat lower energies, since the condition that the CR 
gyroradius becomes comparable to L c in a magnetic field of strength B rms implies 
E ss EJN^ 2 . 



7. Discussion 

Let us now briefly comment on the possible impact of these results for the observation 
of UHECRs. Random magnetic fields are present in the Galaxy, with few /iG strength 
and maximum turbulence scale L max ~ 100 pc, which implies a coherence length 
L c somewhere between 20 and 50 pc, depending on its spectral properties. Random 
fields may also be present on supercluster scales, with much larger coherence lengths, 
L c ~ Mpc, and with strength ~ 10~ 8 -^10~ 7 G. In the Galactic case, the typical energy 
at which large magnification effects can be present, corresponding to that for which 
5 ~ L c /L, is just 

3/2 



£ e *Z41EeV% -M (7-1) 
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Figure 15: Images of the source obtained from the ray-shooting simulation corresponding 
to the previous Figure. Only those with amplification around and above 1/10 are shown in 
the top panel, while in the bottom panels each dot represents 1/100 of the unlensed flux. 



This is larger than the typical energy at which lensing effects associated to the regular 
galactic magnetic field would appear (E/Z ~ 10 EeV). Furthermore, if the random 
fields are concentrated near the galactic plane rather than having a spherical halo 
like distribution, CRs arriving at small galactic latitudes will have associated larger 
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values of L, and hence will suffer lensing effects at higher energies than those arriving 
from higher latitudes. The most important signatures of this would be the likely 
appearance of significant magnification peaks (A int > few) at energies close to E c 
(typically E c /3 < Ej < 2E C ). These peaks will be associated with the appearance 
of the first image pairs, which are very magnified and appear displaced from the 
principal image by an angle ~ L c /L (but the two new images appear on the same 
spot in the sky). This will hence clearly lead to an enhanced signal in a narrow energy 
bin (the typical width of the magnification peak, i.e. ~ 10% of Ef), and can then 
be a source of clustering of events. This is similar than what was previously noticed 
in relation to the regular fields fll2 |, but manifests at somewhat higher energies and 



smaller angular scales (see Eq. ( |2.9| )). In the next Section we actually look for the 
presence of these kind of signatures in the AGASA data above 40 EeV, and find some 
significant hints that could plausibly have their origin in a lensing phenomenon. 

For decreasing energies the number of images increases exponentially, but the 
lensing effects of each one is suppressed because the associated peaks in the spectrum 
become quite narrow. In this way one arrives to the scintillating regime, with many 
images of each source and an overall magnification of order unity. The angular 
extent of this blurred image of the source is given by S rms . The presence of a regular 
component in the magnetic field would further suppress the peaks produced by the 
random component. On the other hand, when for decreasing energies one reaches the 
regime where strong lensing effects associated to the regular field itself are produced, 
these may be somewhat smoothed by the finite extension of the large number of 
subimages just mentioned. In any case, one may understand the effects of the large 
scale regular fields as resulting from the folds produced on a sky which has already 
been corrugated by the action of the random field on a much smaller scale. 

Regarding the extragalactic random fields, the strong lensing effects appear at 
energies 

^Z^eV^J^f 2 ^. (, 2 ) 

which are typically much higher than those associated to the Galactic fields. These 
would hence be relevant for supra GZK energies (E > 10 20 eV) even for CR protons. 

It should also be noticed that galactic scale random fields present around extra- 
galactic sources may be relevant in producing an energy dependent beaming of the 
fluxes which could also lead to interesting features in the observed spectra. 

We certainly look forward to the increased CR statistics that will be accumulated 
at ultra high energies in the near future, and which will allow these effects to be better 
scrutinized. 

8. Epilogue: hints of lensing in the AGASA data? 

An interesting signature of the lensing effects, associated to the high magnification 
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peaks produced when new image pairs appear at the caustic crossings, is the predic- 
tion of an excess of angular clustering of events with similar energies (actually with 
similar rigidities). As we already noticed in fll2"fl , this was strikingly apparent in the 
list of doublets and triplets obtained from the combined data of AGASA, Haverah 
Park, Volcano Ranch and Yakutsk (eight doublets and two triplets within 4° angular 
separation) [|I5|]. Indeed, most of the doublets found there are consistent with the 
two events having the same energy, in one of the triplets the energies are in the ratio 
1:2:4 (as would result e.g. from one event being a proton, the other a He nucleus 
and the other a Be nucleus, all with the same rigidity) while in the second triplet 
two events have similar energies. 

Regarding the analysis of the AGASA data alone [16, 17 1 , involving five doublets 
and one triplet within 2.5° angular separation and E > 4 x 10 19 eV, two of the 
doublets are consistent with having the same energy while two events in a triplet 
also have similar energies. To make these statements more quantitative we have 
studied the correlations among the AGASA published events both in angle and in 
energy, and compared them with simulated events to see the significance of any excess 
observed. 



Following the studies in Jl8|, [17|, [19|, in which one and two dimensional angular 
correlations were analyzed to put in evidence the excess clustering on small angular 
scales and any possible coherent deflection produced by a regular component of the 
galactic magnetic field, we analysed the correlations of the angular separation of the 
events (a = \9 2 — #i|) and their ratio of energies (R = max(E 2 , Ei)/min(E2, -Ei)), 
where Q 12 and E 1)2 are the angular positions and energies of all possible event pairs 
chosen among the AGASA (or simulated) data set. 

Taking several bins in a and R, we defined the density of pairs in the bins 
/(Aa, AR) and plotted the difference between the numbers obtained from the data 
and the corresponding averages obtained from a large set of simulated data and 
normalizing this to the dispersion in the simulated data, i.e. 

Each simulation had the same number of events as the real data within the angular 
cuts performed. Since the exposure of AGASA is uniform in right ascension, no 
cuts were imposed in this angular variable, but to reduce the sensitivity of the anal- 
ysis to the unspecified variation of the exposure with declination (which decreases 
for increasing departures from the latitude of Akeno, which is 35°46'), we only con- 
sidered events with declination in the range [0°,70°], i.e. essentially within ±35° 
of the location of the experiment. This leaves 51 out of the 58 published events 
above 4 x 10 19 eV. We then performed several thousands simulations of sets of 51 
events distributed randomly within the same angular cuts, with an energy spectrum 
dN/dE (x E~ 2J and energies above 40 EeV. 
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Figure 16: Significance of the correlations in angle and energy ratio in the observed 
AGASA data. The right panel corresponds to the first angular bin of the left panel (a < 
2.5°). 



The results are shown in the left panel of Figure [16]. A significant excess of pairs 
is found both at small angles and at similar energies. The right panel shows the one 
dimensional correlation in R for pairs separated by less than 2.5°, which corresponds 
to the first angular bin in the left panel. We see here a more than three sigma excess 
for events with E% ~ E 2 , which reflects the fact that two doublets and a pair in 
the triplet have this property. The other quite significant peaks around R = 3 and 
4 reflect the fact that there are doublets with E 2 /Ei ~ 4.2,3.7 and 3.1, something 
unlikely in a steeply falling spectrum like the one observed and with the present 
statistics. Magnetic lensing of sources with varied composition is thus a plausible 
cause of these peaks. 

Another feature which can be noticed in the observed multiplets is the fact that 
all of them have one or two events with energy near or below 50 EeV, what could be 
indicative of the presence of a threshold energy associated to the formation of caustics 
(with the higher energy events in the pairs corresponding to heavier nuclei). This 
fact also explains the excess observed at similar energies but large angles (15° -^25°) 
in the two dimensional correlation (left panel), which seems to be the result of cross 
correlation between different doublets having similar energies. 
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